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-«.^ , Abstract 

^ ■ 
Q,. We consider a scenario where the strong coupling constant was 

' I changing in the early universe. We attribute this change to a varia- 

;_H ' tion in the colour charge within a Bekenstein-like model. Treating the 

c/3 . vacuum gluon condensate < G^ > as a free parameter , we could gen- 

erate inflation, similar to a chaotic inflation scenario, with the required 
properties to solve the fluctuation and other standard cosmology prob- 
^ ■ lems. A possible approach to end the inflation is suggested. 
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1 Introduction 

The hypothesis that the universe underwent a period of exponential expan- 
sion at very early times is the most popular theory of the early universe. The 

*fellow of CONICET 
tmember of CONICET 



"old" jlj and "new" j21 El inflationary universe models were able to solve the 
"horizon", "flatness" and "structure formation" problems creating, in turn, 
their own problems, some of which the "chaotic" model [ij and its exten- 
sions [Sj could tackle. On the other hand, these kind of cosmological models 
are favored by recent measurements from the Cosmic Microwave Background 
(CMB) jSmi and the WMAP results |8j. Usually, inflationary models are 
all based on the use of new fundamental scalar flelds, the "inflatons" , which 
can not be the Higgs fields of ordinary gauge theories. Later, other possible 
alternatives to inflationary cosmology were proposed. Rather than chang- 
ing the matter content of the universe, new concepts were adopted such as 
non-commutative space-time quantizations [9^, brane settings of a cycling 
universe PU] , and a change in the speed of light in the early universe [TT| IT^ 
or, more generally, a varying flne structure constant ttem [IH]- In [HI; we 
have generalized Bekenstein model ^Hj for the time variation of aem to QCD 
strong coupling constant as and found that experimental constraints going 
backward till quasar formation times rule out as variability. In this letter, we 
discuss how our "minimal" Bekenstein-like model for as, when implemented 
in the very early universe, can provide a realization of inflation driven by 
the trace anomaly of QCD energy momentum tensor. Inflationary models 
driven by the trace anomaly of conformally coupled matter flelds are treated 
in the literature |16j while in our toy model we concentrate on the gauge 
flelds of QCD. Assuming the universe is radiation-dominated at early times 
and that vacuum expectation values predominate over matter densities, we 
flnd, with suitable values of the gluon condensate < G^ > far larger than 
its present value and the Bekenstein length scale / smaller than the Planck- 
Wheeler length scale Lp, that our model for inflation is self-consistent with 
acceptable numerical results to solve the fluctuation and other problems. 
We flnd that the QCD lagrangian with a changing strong coupling constant 
leads naturally to a monomial quadratic potential like the chaotic scenario. 
However, while the large values of the inflaton matter fleld plague the latter 
scenario, they just amount in our model to a reduction of the strong charge 
by around 10 times during the inflation. We shall not dwell on the possible 
mechanisms by which the gluon condensate could have decreased to reach 
its present value, and which might be necessary in order not to have exotic 
relics. Rather we wish to concentrate on the conditions we should impose 
on our model to be interesting. We hope this phenomenological approach 
could prompt further work on a possible connection between time-varying 
fundamental "constants" and inflationary models. 



2 Analysis 

We follow the notations of ^^. Therein, we used the QCD Lagrangian with 
a varying coupling "constant" 



^QCD 






2/2 e2 2 ^ 

where the metric used for raising and lowering indices is the R-W metric, / is 
the Bekenstein scale length, e{x) is a scalar gauge-invariant and dimensionless 
field representing the variation of the strong coupling "constant" g{x) = 
goe{x), D^ = d^ — igoe{x)Af^ is the covariant derivative and where the gluon 
tensor field is given by: 

G;. = -^meAt)-dAeA;)+goe'r'^A'^At] (2) 

We assume homogeneity and isotropy for an expanding universe and so con- 
sider only temporal variations for 0:5 = ^^^ = asoe^{t). We assume also, 
rather plausibly in the radiation-dominated early universe, negligence of mat- 
ter contribution to get the following equations of motion 

i^\,- gor''Gt'Al + J2go4^t^Y^ = (3) 

^ / 



(a^i)- = ^{G^) (4) 

where a{t) is the expansion scale factor in the R-W metric. 

Computing the canonical energy-momentum tensor QiQj,.\ d^4'i ~ g^^L 
we get a non gauge invariant quantity. This may be cured by a standard 
technique |17j amounting to a subtraction of a total derivative and hence 
not changing the equations of motion. We subtract the total derivative 
AT°^ _ Q^ f — a-^j[a.p\ |.Q gg|.^ with the use of equation (jHl), the gauge- 
invariant energy momentum tensor 



j.a(3 ^ GTGf + ^5:^(^)7(-D^V^^^ Id^ed^e 
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\g':G% + Y.^'^f\ii^D, - Mf)i;^f^ - i^^ 



(5) 



Since we assume radiation dominated era during the very early universe 
we can concentrate on the gauge fields and neglect the matter fields contri- 
bution and so we decompose our energy momentum tensor into two parts: 
the gauge part and the e-scalar field part 

T^p = n, + T^, (6) 



^au^R 9aP 



A CL ^ fJiU 



1 dgedne _ 



:Ml 



2/2 e 



2 



Here all the operators are supposed to be renormalized and it is essential in 
the infiationary paradigm that quantum effects are small in order to get small 
fluctuations in the CMB, which will be seen to be the case in the model. 

The contribution of the scalar fleld to the energy density p^ = Tqq and to 
the pressure T^^ = g^j Pe can be computed and we get 

Pe =-M'f= V. (7) 

On the other hand, the gauge field contribution Tf^ can be decomposed 
into trace and traceless parts. In such a way, we can write the corresponding 
energy mass density pg = Tqq and the pressure Pg, like in "ordinary" QCD, 
as a sum of two terms: 



P9 = pI + P^ (8) 

Pa = pI+P^ (9) 

where p^ is the density corresponding to the "traceless" part of the gauge 
field satisfying 

Pi = ^Pl (10) 

while the trace part of the gauge field energy-momentum tensor is propor- 
tional to gai3 and behaves like a 'cosmo logical constant' term. Thus, the 
corresponding equation of state reads: 

Pj = -Pj (11) 

We shall assume, here, that the trace anomaly relation for T^^ is the same 
as the "ordinary" QCD trace anomaly. This can be justified because the 
energy-momentum tensor T^o is identical in form to "ordinary" QCD and 
that the trace anomaly which involves only gauge invariant quantities should, 
by dimensional analysis, be proportional to G^. However, we have checked 



that a change in the numerical value of the proportionality factor would not 
alter our conclusions. Thus we take, up to leading order in the time-varying 
coupling "constant" as = a^o^^; ^^^ relation [T8] : 



^r = p. - 3p. = -^GTc;:, (12) 

Again, neglecting matter contribution during the radiation dominated era, we 
limit our gluon operator G'^^G^^,^ matrix elements to its condensate vacuum 
expectation value < G^ >, and so we get: 



Hence, we can write the total energy mass density as 

P = P^ + Pj + P. (14) 

and equation (fTT|) would suggest, in analogy to ordinary inflationary models, 
that the QCD trace anomaly could generate the inflation. For this, let us 
assume that the "trace-anomaly" energy mass density contribution is much 
larger than the other densities: 

p'^»Pe,pl => P^P^ (15) 

Then, equation (fTSj) tells that the vacuum gluon condensate < G^ > should 
have a negative value which is not unreasonable since the inflationary vacuum 
has "strange" properties. In ordinary inflationary models, it is filled with 
repulsive-gravity matter turning gravity on its head pH]- This reversal of 
the vacuum properties is reflected, in our model, by a reversal of sign for the 
vacuum gluon condensate. 

Now we seek a consistent inflationary solution to the FRW equations in 
a flat space-time: 

d 2 SttGn , . 

(-) = ^^— P (16) 

- = 5— P + 3p 17 

a 6 

where Gn is Newton's constant. The first FRW equation with (|T!?jl will give 
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On the other hand, the equation of motion (eqEj) of the scalar field can be 
expressed in the following way: 

- + m--{-f = (19) 

This equation differs from the ordinary 'matter' inflationary scenarios in the 
term (-)^. However, for "slow roll" solutions we neglect the terms involving 
f and {\f to get 

2>He =!^^e = -V'{e) (20) 

which is the same as the "slow roll" equation of motion of the inflaton in 
ordinary scenarios. In our model, the "slow roll" condition can be written 
as: 

5=1-^1 =^(^)H «1 (21) 

This gives us the first hint that Bekenstein hypothesis (Lp < /) might not 
survive. Now, we set e/, the value of e-field at the end of inflation t/, to 1 
so that the time evolution of the strong coupling terminates with the end of 
inflation. We expect for "slow roll" solutions that ej, the value of e at the 
initial time of inflation tj corresponding to when the CMB modes freezed out, 
to be of order 1. Taking this into account, and in order that the changes of 
the Hubble constant and the energy mass density are not very large during 
the inflation, we assume the gluon condensate value < G^ > to stay approxi- 
mately constant during much of the inflation. Then we get a linear evolution 
in time for the e- scalar field 

e{t) = e.-:T^7j7^-^TJ^iU'GT\<G'>\'/'{t-U) (22) 

3'^/^(aso)^ Lp 

and we have, as in chaotic scenarios, a simple quadratic potential: 

V{e) = LK^^2 (23) 

One can make explicit the correspondence between our model with e- 
scalar field and the chaotic scenario with an 0-inflaton matter field. Look- 
ing at equation ()13|) and comparing with p = V{<f>) in ordinary inflationary 
models, we see that the gluon condensate plays a role of a potential for the 
"inflaton" -e field. Comparing equations ()18p and l^20\i with the corresponding 
relations in ordinary inflationary models: 

H^ = ^««V(^) (24) 

3H^ = -V'(0) (25) 



and remembering that the dimension of the matter scalar field (j) is one, we 
^ and oc |. 



have V oc ^ and </> oc y, and we can find the relations between (e, V{e)) and 



(0,V(<^)): 

0=^e with y=^ (26) 

lv{e) =V(0)= ^''^f ^' 0^ (27) 

3 Results and Conclusion 

Now, we check that our model is able to fix the usual problems of the standard 
(big bang) cosmology. First, in order to solve the "horizon" and "fiatness" 
problems we need an infiation -jA- of order 10^^ implying an infiation period 
At = t f — ti such that 

H At r^ 65 (28) 

Furthermore, it should satisfy the constraint 

10-^°s < At < 10-^°s (29) 

so that not to confiict with the explanation of the baryon number and not to 
create too large density fiuctuations 22^ 23j. The bound 10~^°s corresponds 
to the time, after the big bang, when the electroweak symmetry breaking 
took place. Presumably, our infiation should have ended far before this time. 
Thus, from equations (J28|l . (J29|l and |T8|) we get the following bounds on 
I < G^ > |: 

3 X lO'^GeV^ < e| < G^ > |V2 < 3 X 10^'^GeV^ (30) 

In order to determine e,, we have ^^^ = ^ ~ —^h) — ~ m'^v' which gives, 
using equations (|T!?jl and pUjl . the relation: 

65 ~ln^= (^)^^(e?-l) (31) 

Next, comes the "formation of structure" problem and we require the frac- 
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tional density fiuctuations at the end of infiation to be of the order 
10^^ so that quantum fiuctuations in the de Sitter phase of the infiationary 
universe form the source of perturbations providing the seeds for galaxy for- 
mation and in order to agree with the CMB anisotropy limits. Within the 
relativistic theory of cosmological perturbations |24j, the above fractional 



density fluctuations represent (to linear order) a gauge-invariant quantity 
and satisfy the equation 

p 

where 5M represent the mass perturbations. 

The initial fluctuations are generated quantum mechanically and, at the 
linearized level, the equations describing both gravitational and matter per- 
turbations can be quantized in a consistent way J2S|. The time dependence 
of the mass is reflected in the nontrivial form of the solutions to the mode 
equations and one can compute the expectation value of field operators such 
as the power spectrum and get the following result for the initial mass per- 
turbation 121123 

whence 

In order to evaluate (p + p)\ti we use the energy conservation equation: 

p + 3{p + p)- = (35) 

a 

and after substituting p r^ p^ we get 

In fact, the energy conservation equation can be used to solve for p^ and we 
could check that 

PliPl)-Pe{Pe) - 6xp^{6xpT) (37) 

where 5 = \-^\ ~ ^(r~)^ ^^^ ^'-'' ^^^^ the "slow roll" condition ()21|1 is 
satisfied, our solution assuming the predominance of the "trace-anomaly" 
energy mass density is self-consistent. Substituting equation (jHH|l in (jH^ and 
using equation (plj) we get 

^ ~ 9y^a5olOV|<G^>|^G^ (38) 



Hence, taking G^ ~ 10 ^^GeV ^ we obtain 

/ |<G'2>|5 



rvj 



e' (39) 



Lp 1034^61/2 

and combining this last result with (jHUj) . we get 

10"^^ « -— < 10^ (40) 

e Lp 

The "slow roll" condition (j'^lj) is consistent with the upper bound, while the 
lower bound restricts e^ not to be too large. 

On the other hand, it is possible to calculate the spectral index of the 
primordial power spectrum for a quadratic potential as follows: 



n - 


- 1 = -4?7 


where rj = 


Mil y 2 
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Lp 


'h^ 


(41) 


and we find: 




n = 1- 








(42) 


The inflation would end 


(e/ = 1) when the "slow roll 


" parameter r] = 


6 = 1. 


We should evaluate the 


QCD coupling constant as^ 


(/^) = 


47r 

/3oln(-^) 
QCD 


■ at an 



energy scale corresponding to the inflationary period. We take this to be 
around the GUT scale ~ lO^^Ge^ and /?o = ll — |?^/ = 7 whereas the weak 
logarithmic dependence would assure the same order of magnitude for as^ 
calculated at other larger scales. With Aqcd ~ 0.2GeV fIU\ we estimate 
aso ~ 0.025, and so we get 

i^? ~ 10-^ (43) 

Lp 

This is in disagreement with Bekenstein assumption that Lp is the shortest 
length scale in any physical theory. However, it should be noted that Beck- 
enstein's framework is very similar to the dilatonic sector of string theory 
and it has been pointed out in the context of string theories 26j that there 
is no need for a universal relation between the Planck and the string scale. 
Furthermore, determining the order of magnitude of -^ is interesting in the 
context of these theories. 

From (jHH), we have e^ ~ 11, and then using ()42|1 . we have n = 0.97 which 
is within the range of WMAP results ^ . 

9 



Thus, we see that our model reproduces the results of the chaotic infla- 
tionary scenario. However, the shape of the potential was not put by hand, 
rather a gauge theory with a changing coupling constant led naturally to it. 
Moreover, in typical chaotic models, the infiaton field starts from very large 
values {(pi ~ 15Mpi) and ends at around IMpi. One might suspect whether 
field theory is reliable at such high energies. Nonetheless, this problem is 
absent in our model since the large values have another meaning in that they 
just refer to a reduction of the strong coupling by around 10 times during 
the inflation. 

Furthermore, chaotic inflations get a typical reheating of order Trh ~ 
lO^^GeV, and one might need to worry about the relic problem. Similarly, 
equation (jHIJjl leads in our model to a gluon condensate | < G^ > |j ~ 
lO^^GeV^ at the start of inflation. From equation (j221); we see that this 
corresponds to an inflation time interval At ~ 10~^^s satisfying the con- 
straint (j2ni)- If the gluon condensate stays constant, as we assumed in our 
analysis, we will have the same reheating temperature as in chaotic models 
(Trh ~ p{tfY^'^). However, we should compare this value for < G^ > with its 
present value renormalized at GUT scale ~ lO^^GeV which can be calculated 
knowing its value at IGeV ^1] and that the anomalous dimension of asG'^ 
is identical ly zero. We get 

< G\now, li ~ IQ^^GeV) >~ IGeV^ (44) 

which represents a decrease of 62 orders of magnitude. 

This can give us the following possible picture for an exit scenario. Lack- 
ing a clear theory for the non-perturbative dynamical gluon condensate, we 
consider its value | < G^ > | depending on energy, and thus implicitely on 
cosmological time, as given by the standard RGE which turns it off logarith- 
mically at high energy. However, we can furthermore assume the condensate 
value to depend explicitely on time during inflation: 

<G\E,t)>=<Gl{E{t))>f{t) 

where < Gq{E) > is the piece determined by the RGE. In order that our 
model be consistent, the value of < G^ > at the start of inflation should be 
very huge and negative. The unknown function f{t) should be such that it 
varies slowly during most of the inflationary era, to conform with an approx- 
imately constant value of < G^ >, while at the end of inflation it causes a 
drastic drop of the condensate value < G^ > to around zero. The energy re- 
lease of this helps in reheating the universe, while reaching the value leads 
to a minute "trace- anomaly" energy mass density (equation II H|l ending, thus, 
the inflation. The other types of energy density would contribute to give the 

10 



gluon condensate its 'small' positive value of (jii|) . and the subsequent evo- 
lution is just the standard one given by RGE. Surely, this phenomenological 
description needs to be tested and expanded into a theory where the concept 
of symmetry breaking of such a phase transition for the condensate < G^ > 
provides the physical basis for ending the inflation. Nonetheless, with a test 
function of the form f(t) = — /3^tanh^(e — 1) with /3 ~ 10^^, one can inte- 
grate analytically the equation of motion, and in "slow roll" regime we have 
e = 1 + Arcsinh{exp[—a(3t]) with a(] ~ lO^^GeV. The graph in Fig. Qshows 
the time evolutions of the condensate and the e-field, which agree with the 
required features. This example is meant to be just an existence example, 
and the temporal dependence of the condensate could be of complete differ- 
ent shape while the whole picture is still self-consistent. The issue demands 
a detailed study for the condensate within an underlying theory and we do 
not further it here. We hope this work will stimulate interest in the subject. 




time 
t i t f 



Figure 1: Temporal evolution of the condensate | < G^ > | (thick line) and 
the e-field (thin line), for the choice /(t) = -f]'^tanh^{e - 1). The < G^ > 
scale has been adapted so that to visualize both graphs together. 
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